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DATA  MANIPULATION  PACKAGE  USING  PV-WAVE 


INTRODUCTION 

The  Data  Manipulation  Package  (DMP)  is  a  computer  program  that  was  written  in  the 
PV~WAVE  [1]  command  language.  PV-WAVE  is  a  high  level  data  analysis  and  graphics 
programming  language  that  is  useful  in  data  manipulation  and  graphical  presentations.  Unfortunately,  to 
use  the  PV-WAVE  language  effectively,  the  user  must  be  familiar  with  virtually  hundreds  of 
commands,  functions,  procedures  and  system  variables.  The  DMP  program  was  developed  to  facilitate 
the  use  of  PV-WAVE  by  non-programmers  by  providing  a  menu-driven  user  interface.  The  capabilities 
supported  by  DMP  are  subdivided  into  four  categories:  plotting,  curvefitting,  data  manipulation,  and 
statistical  calculations.  Each  of  these  categories  has  many  options  that  are  described  in  detail  in  this 
report 

DMP  is  written  to  be  as  portable  as  the  PV-WAVE  program  itself.  DMP  was  written  for  use  on 
Digital  Equipment  Corporation’s  VMS  machines,  as  well  as  most  UNIX  based  workstations. 

This  report  is  not  a  user’s  manual,  but  rather  serves  as  a  general  overview  for  the  DMP  program. 
The  first  section  outlines  the  capabilities  of  DMP  and  the  second  section  provides  general  operating 
instructions.  The  appendices  contain  information  regarding  the  organization  of  DMP  along  with  the 
subroutine  descriptions.  The  computer  code  is  available  to  all  interested  persons  by  requesting  a  copy 
from  the  following  address: 

Naval  Research  Laboratory 
Underwater  Sound  Reference  Detachment 
Attn:  Code  5916  (C.  M.  Siders) 

P.O.  Box  568337 
Orlando,  FL  32856-8337 

PROGRAM  CAPABILITIES 

DMP  addresses  the  data  analysis  and  graphical  presentation  requirements  deemed  most 
necessary  in  the  everyday  operation  of  a  scientific  laboratory  using  the  PV-WAVE  command  language. 
PV-WAVE  can  perform  many  operations  that  are  not  included  in  the  DMP  interface;  however,  these 
options  could  be  added  to  DMP  if  necessary.  Figure  1  is  a  diagram  illustrating  the  basic  subdivisions  of 
DMP;  it  does  not  indicate  the  program  flow  but  organizes  the  capabilities  of  DMP  in  a  logical  manner. 
The  first  entry  in  the  diagram  in  Fig.  1  is  DATA  ENTRY.  There  are  three  ways  to  bring  data  into  DMP: 
read  data  from  a  file,  generate  data  within  the  DMP  program,  and  digitize  data  from  a  compatible  tablet 
Data  entry,  of  course,  must  be  performed  before  any  other  function.  Once  the  data  is  entered  it  then 
must  go  through  an  "identification"  phase.  To  identify  the  data  simply  means  to  select  the  column  or 
columns  necessary  for  the  desired  operation.  For  example,  if  a  cartesian  plot  is  required,  the  user  must 
select  the  two  columns  that  will  be  identified  as  the  abscissa  and  the  ordinate. 
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Fig.  1  -  Diagram  depicting  DMP  capabilities 

After  the  data  has  been  loaded  and  selected,  the  user  has  several  options  as  shown  in  the 
"ANALYSES"  entry  in  Fig.  1.  The  output  for  these  procedures  can  be  a  visual  display  or  a  newly 
created  file,  or  both.  A  description  of  the  functionality  of  each  of  the  specific  routines  will  be  described 
in  more  detail  below. 
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Data  Entry 

One  way  of  making  data  available  for  the  DMP  program  is  to  load  data  from  a  file.  This  file 
should  be  columns  of  numbers  (data  set)  stored  in  ASCII  format.  The  numbers  can  be  separated  by 
commas  or  any  number  of  blanks.  Each  row  of  data  in  the  file  must  have  the  same  number  of  columns. 
The  program  queries  the  user  for  the  name  of  a  file  to  be  read.  If  the  user  is  unsure  of  the  filename(s)  to 
be  loaded,  he  may  press  the  ’L’  key  to  list  all  files  in  the  current  directory  with  a  .dat  extension.  Once 
the  filename  has  been  entered.  DMP  loads  the  file  and  prints  the  number  of  rows  and  columns  in  the  file 
on  the  screen.  The  user  then  has  the  opportunity  to  load  more  data  or  return  to  the  main  menu. 

Data  can  also  be  generated  with  the  DMP  program  illustrated  by  the  "GENERATED"  block  in 
Fig.  1.  This  routine  queries  for  a  beginning  value,  an  ending  value,  and  the  number  of  points.  A  single 
column  of  data  is  generated  and  can  be  stored  permanently  in  a  file.  This  data  generation  option  may 
seem  limited;  however,  with  the  available  data  manipulation  capability,  the  user  can  execute  a  variety  of 
operational  forms.  These  general  forms  will  be  listed  in  the  "DATA  MANIPULATION"  section  in  a 
latter  part  of  this  report. 

The  final  method  to  input  data  is  by  digitizing  points  from  a  compatible  tablet.  DMP  requires 
the  user  to  mark  the  lower  left-hand  and  upper  right-hand  comers  of  the  drawing  area.  The  lower  left 
point  should  be  placed  at  the  minimum  value  of  the  abscissa  and  ordinate.  The  upper  right  point  should 
be  placed  at  he  maximum  value  of  both  the  abscissa  and  ordinate. 

Data  Selection 

Once  all  the  data  files  are  loaded,  the  user  is  asked  to  "IDENTIFY"  the  data.  To  identify  the 
data  means  to  select  the  filename  and  column(s)  that  will  be  used  for  further  procedures.  All  the 
PROCEDURES  shown  in  Fig.  I  require  that  either  one  or  two  columns  of  data  be  selected. 

Analyses 

There  are  four  main  sections  that  will  be  discussed  in  detail:  Plotting,  Curvefitting,  Data 
Manipulation  and  Statistical  Evaluations.  When  a  new  data  set  is  created  from  one  of  these  sections,  this 
does  not  imply  that  the  data  is  permanently  saved  to  a  file.  The  name  of  this  data  set  is  used  in  the 
current  session  only,  unless  the  user  specifically  chooses  to  save  data  to  a  file.  If  a  graphic  display  is 
available,  it  may  be  viewed  on  the  CRT  screen  or  sent  to  a  hard  copy  device.  Appendix  A  details  how 
the  graphic  output  routines  can  be  customized  to  the  user’s  needs.  Also,  data  can  go  from  one  routine  to 
another.  For  example,  loaded  data  can  be  curve  fitted,  then  manipulated  and  finally  plotted. 

Plotting 


The  plotting  capability  in  DMP  was  developed  to  give  the  amount  of  flexibility  needed  for 
publication  quality  displays  without  enduring  extensive  interactive  computer  sessions.  This  is  done  in 
DMP  with  menu-driven  screens  in  which  all  the  options  have  intrinsic  defaults  which  allow  the  user  to 
produce  plots  quickly  by  keeping  the  defaults;  however,  DMP  allows  the  user  to  customize  wherever 
needed.  The  customization  parameters  are  used  throughout  the  entire  DMP  session.  There  are  two  basic 
plotting  capabilities:  Polar  Plotting  and  Cartesian  Plotting. 

Cartesian  Plotting.  Three  forms  of  cartesian  plotting  include  single  plot,  overplot,  and  multiple 
plot.  "Single  plot"  produces  one  curve  on  a  set  of  axes;  "overplot"  produces  several  curves  on  the  same 
axes;  and  "multiple  plot"  produces  two  to  four  or  smaller  single-curved  plots  on  a  page. 
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)  Line/ symbol  type:  Solid 
)  Line  color:  White 

)  output  device:  VT340  Screen 

)  Printer  parameters  menu 
)  Mouse  Utilities 


R . )  Return  to  previous  menu 
Q.)  Quit  from  plotting  package 


Enter  choice: 


Fig.  2  -  Standard  plot  menu  viewed  in  an  interactive  DMP  session. 

Single  Plot.  Figure  2  shows  the  many  customizing  features  available  in  the  Single  Plot  option. 
Any  text  to  the  right  of  the  colon  is  a  default  value.  If  a  user  simply  wants  a  quick  plot,  he  chooses 
entry  l)Plot.  An  autoscaled  plot  with  no  labels  would  appear  on  the  screen.  However,  if  customization 
is  necessary,  options  2  through  12  can  be  selected.  Options  4  and  7  allow  the  user  to  choose  between 
linear  and  logarithmic  scaling.  Options  5  and  8  allow  the  user  to  choose  between  autoscaling  the  data  or 
manually  entering  minimum  and  maximum  values  for  the  abscissa  and  ordinate.  Only  the  data  in  the 
prescribed  range  will  be  plotted.  The  line  type  (Option  9)  can  either  be  a  line,  a  symbol,  or  a 
combination  of  the  two.  Option  12  allows  the  user  to  change  plot  dimensions  and  orientation  for  output 
sent  to  printers(landscape  or  portrait).  Portrait  mode  is  the  default  in  which  the  top  of  the  page  is  8  1/2 " 
wide.  Mouse  utilities  (option  13)  requires  that  a  mouse  and  mouse  driver  be  installed.  This  option 
provides  a  cross-hair  that  moves  and  executes  under  control  of  the  mouse.  When  the  cross  hairs  are 
placed  close  to  a  point  on  the  curve  and  the  mouse  is  clicked,  the  values  of  the  ordered  pair  will  be 
returned.  The  user  is  also  able  to  zoom  on  a  portion  of  the  plotting  area  and  detemrine  local  minimums 
and  maximums  within  the  restricted  zone.  Auxiliary  instructions  are  given  on  the  last  line  of  the  screen. 
An  example  of  a  single  plot  is  given  in  Fig.  3. 
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Test  #1 


Fig.  3  -  Example  of  a  "Single  Plot"  option. 

Overplot.  The  overplot  option  requires  that  all  the  curves  to  appear  on  the  same  plot  be  marked 
by  a  toggle  feature  in  DMP.  The  only  other  additional  feature  with  this  option  is  that  the  user  can  locate 
a  legend  anywhere  in  the  plotting  area  by  entering  the  coordinates  of  the  lower  left  hand  comer.  The 
legend  consists  of  the  line  type  along  with  a  short  description  as  shown  in  Fig.  4.  In  order  to  access  the 
legend,  the  user  must  enter  the  "Line  Definition  Menu".  An  example  of  an  overplot  with  legend  is 
shown  in  Fig.  5. 
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Fig.  4  -  DMP  menu  illustrating  legend  and  line  characteristics. 
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Fig.  5  -  Example  of  "Overplot"  with  legend. 

Multiple  Plot.  The  multiple  plot  option  allows  the  user  to  place  several  (2-4)  smaller  versions  of 
a  single  plot  on  a  page.  These  plots  are  selected  by  a  toggle  feature  in  DMP  similar  to  that  in  the 
overplot  option.  The  user  can  select  the  portrait  or  landscape  mode  and  can  choose  a  row  major  and 
column  major  priority.  These  priorities  allow  the  user  to  arrange  the  plots  by  row  first  or  by  column 
first.  Multiple  plots  in  row  major  priority  are  shown  in  Fig.  6.  Note  the  number  in  the  upper  left-hand 
comer  shows  the  order  in  which  the  plots  were  selected. 
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Fig.  6  -  Example  of  "Multiplot  Option”  in  row  major  priority. 

Polar  Plotting.  The  "Polar  Plotting"  option  provides  the  capability  of  overplotting  as  shown  in 
Fig.  7.  The  abscissa  and  ordinate  for  a  polar  p!  jt  must  be  the  angle  (default  is  radians)  and  the  radius, 
respectively.  PMP  queries  the  user  for  the  number  of  rings  and  the  radial  difference  between  rings.  The 
user  also  has  the  option  of  setting  any  negative  ndial  values  to  zero  in  the  data  file. 
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Fig.  7  -  Example  of  polar  plot  implementing  overplot  feature 


Curvefitting 


The  user  has  a  variety  of  functional  forms  to  choose  from  for  data  curvefitting.  All  these  forms 
allow  the  user  to  generate  data  beyond  the  range  of  the  data  that  is  being  curvefitted.  A  Goodness  of  Fit 
Criterion  is  calculated  and  displayed  for  all  the  fitted  data  except  the  spline  fit.  The  value  for  this 
parameter  ranges  from  zero  to  one.  The  closer  the  parameter  is  to  one  the  better  the  fit  The  curvefitting 
forms  include  polynomial,  rational  fraction,  exponential,  logarithmic,  and  cubic  spline.  The  data  used 
for  the  rational  fraction,  exponential,  and  logarithmic  curvefitting  routines  is  automatically  scaled  to 
avoid  any  numerical  overflow  problems. 

The  polynomial  fit  includes  all  terms  up  to  the  eighth  order.  The  user  may  choose  an  order  less 
than  eight,  however,  for  each  order  all  preceding  terms  exist.  The  rational  fraction  fit  has  the  form: 


C.  +  C.x  +  C,x^  +  C.x^  +  C_x4  +  C.x^  +  C_x^ 

1 _ 2 _ 3 _ 4 _ 5 _ 6 _ 7 

C8  *  V  *  C10x2  *  Cll*3  *  C12x4  +  C13x5  *  C14x6 


(1) 


Unlike  with  the  polynomial  fit,  the  user  is  allowed  to  omit  any  terms.  The  exponential  and 
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logarithmic  fit  have  the  following  torms  respectively: 

c  *  exp(-bx)  +  d  (2) 

and 

a  *  ln(  bx)  +  d.  (3) 

The  last  curvefitting  technique  is  a  cubic  spline  interpolation.  This  method  is  useful  in  generating  a 
smooth  curve  through  a  set  of  data  points  without  being  concerned  with  fitting  to  a  particular  form. 
Figure  8  is  an  example  of  a  cubic  spline  fit  through  a  given  set  of  points. 


Fig.  8  •  Example  of  cubic  spline  interpolation  for  a  given  set 


Data  Manipulations 

The  DMP  program  allows  the  user  to  create  a  column  of  data  by  operating  on  one  or  two 
existing  columns  of  data.  The  list  of  symbols  supported  are  shown  in  Fig.  9. 
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SYMBOLS: 

A*F  :  addition  between  two  columns  lof  equal  size) 
A-B  :  subtraction  of  two  columns  tot  equal  size) 

3 *A  :  multiplication  of  a  constant  with  a  column 

A/5  :  division  of  a  column  by  a  constant 

BA2  :  raise  each  element  in  a  column  to  a  power 


Fig.  9  -  A  list  and  description  of  the  symbols  available  in  DMP. 

There  are  advanced  operations  that  can  be  implemented.  Figure  10  is  a  comprehensive  list  of  functions 
that  are  available. 


BUILT  IN  FUNCTIONS: 


ABStX) 

ACOS(x) 

ALOG(x) 

ALOGIO(X) 

ASIN(X) 

AT AN ( X ) 

A  VG  ( X ) 
COSlX) 

COSH (x) 

EXP (X) 
FFTfx, dir) 


MAX (x) 
MEDIAN (X) 
MIN  ( X ) 

SIN  (X) 
SINH(X) 
SQRT(X) 
SUM (X) 

TAN  (X) 
TANH(X) 


<-  compute  the  absolute  value  of  each  element  in  x 
<-  compute  the  arcosine  of  each  element  in  x 
<-  compute  the  natural  logarithm  of  each  element  in  x 
<-  compute  the  base  10  logarithm  of  each  element  in  x 
<-  compute  the  arcsine  of  each  element  in  x 
<-  compute  the  arctangent  of  each  element  in  x 
<-  compute  the  average  of  the  values  of  all  elements  in  x 
<-  compute  the  cosine  of  each  element  in  x 
<-  compute  the  hyperbolic  cosine  of  each  element  in  x 
<-  compute  the  value  eAx  for  each  element  in  x 

<-  compute  the  fast  fourier  transform  of  x.  A  negative  value  for 

dir  transforms  from  the  space  domain  to  the  frequency  domain.  A 

positive  values  for  dir  transforms  from  frequency  to  space 

<-  locates  the  largest  value  in  x 

<-  compute  the  median  value  in  x 

<-  locates  the  smallest  value  in  x 

<-  compute  the  sine  of  each  element  in  x 

<-  compute  the  hyperbolic  sine  of  each  element  in  x 

<-  compute  the  square  root  of  each  element  in  x 

<-  sums  all  values  in  x 

<-  compute  the  tangent  of  each  element  in  x 
<-  compute  the  hyperbolic  tangent  of  each  element  in  x 


Fig.  10  -  A  list  and  description  of  the  functions  available  in  DMP. 


The  "data  manipulation"  option  is  performed  on  a  "selected"  set  of  data,  where  the  first  column 
selected  is  labeled  "a”  and  the  second  column  selected  is  labeled  "b."  The  data  "selection"  process  could 
have  occurred  earlier  in  the  same  DMP  session  under  a  different  procedure.  For  example,  if  the  user 
wants  to  add  3sin  of  the  first  column  of  the  selected  data  to  the  exponential  of  the  second  column  of  the 
selected  data  he  would  type  the  following  command:  3  *  sin  (a)  +  exp  (b). 

In  the  data  manipulation  procedure,  the  user  can  use  either  column  "a"  or  column  "b",  or  both. 
If  only  one  column  is  selected,  then,  of  course,  that  column  would  be  labeled  column  "a."  The  output  of 
this  data  manipulation  is  temporarily  stored  in  a  user  defined  name.  This  data  can  also  be  permanently 
saved  to  a  file. 

Statistical  Parameters 

The  statistical  plotting  package  displays  a  frequency  histogram  (or  relative  frequency  histogram) 
of  data  and  calculates  the  following  parameters:  mean,  median,  mode,  range,  minimum,  maximum,  and 
standard  deviation.  These  parameters  are  printed  to  the  right  of  the  histogram.  The  user  also  has  the 
option  of  displaying  these  parameters  in  graphical  form  by  overlaying  them  on  the  histogram  as  shown 
in  Fig.  1 1.  This  option  requires  that  one  column  of  data  be  selected. 
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Fig.  11  -  An  example  of  a  histogram  with  respective  statistics. 

INTERACTIVE  SCREEN  ORGANIZATION 


Mean  :  46.2160 
Median:  36.0000 


Made  :  5.00000 


Std  Dev:  37.5478 


Range  :  142.000 
Min  :  0.000000 
Max  :  142.000 


This  section  gives  the  user  an  overview  of  the  actual  computer  screen  that  may  be  observed  in  a 
typical  interactive  session.  Figure  12  illustrates  the  hierarchy  of  menu  screens  in  a  block  diagram 
format.  Not  all  of  the  screens  are  included  in  Fig.  12  because  either  they  do  not  contribute  to  the 
understanding  of  DMP  or  they  were  already  described  in  the  text  The  words  in  each  of  the  boxes  in  Fig. 
12  are  the  actual  titles  of  the  interactive  screens  seen  by  the  user. 
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Fig.  12  -  Hierarchy  of  interactive  menu  screens 
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Appendix  A 

OUTPUT  DEVICE  CUSTOMIZATION 


This  appendix  contains  information  useful  in  customizing  DMP  to  support  output  devices  not 
currently  available.  The  devices  supported  by  DMP  include  ReGIS  terminals,  X  windows  terminals, 
TEK4010  terminals,  Talaris  or  QMS  QUIC  printers,  and  Postscript  printers.  Three  PV-WAVE  routines 
must  be  modified  when  adding  a  new  printer  definition  which  includes  change_device.pro, 
setup_device.pro,  and  reset_device.pro. 

Change_device.pro  is  a  routine  that  provides  the  user  with  a  list  of  output  devices  from  which  to 
select.  This  routine  returns  a  number  in  the  variable  "device"  (via  the  devcom  common)  that  contains 
the  user’s  selection. 

The  routine  setup_device.pro  redirects  graphics  to  the  correct  output  device  with  some  minimal 
formatting.  The  formatting  uses  three  parameters  (wx,  wy,  and  orient)  that  are  passed  via  a  common 
statement.  The  variables  wx  and  wy  are  the  length  in  inches  of  the  abscissa  and  ordinate  axis 
respectively.  Orient  determines  if  a  display  is  landscape  or  portrait  mode.  Orient  is  set  equal  to  1  if  the 
landscape  mode  is  desired.  If  orient  is  set  to  any  other  value  the  graph  will  be  in  portrait  mode. 

The  third  routine,  reset_device.pro,  sends  the  file  to  the  proper  queue.  This  routine  also  sets  the 
plot  device  back  to  its  default  value  and  closes  any  files  created. 

If  it  is  necessary  to  add  a  device,  be  sure  to  add  a  unique  device  number.  Users  can  change  the 
order  of  the  devices  by  changing  the  sequence  given  in  the  array  variable  "device_number"  at  the  top  of 
the  program  change_device.pro. 


pro  change_device,  polar_option 


common  devcom,  device,  device_name,  wx,  wy,  orient 

device_names  =  ['VT340  Screen',  'QMS  Printer',  'X  Window',  $ 

’Tek4014  Screen',  'VT340  File',  'QMS  File',  $ 

'Tek4014  File',  ' HPGL  File',  'Pre-printed  format',  $ 
'PostScript  Printer',  'PostScript  File'] 

;  device_number  defines  the  order  that  the  devices  will  appear  in  the  menu. 
;  This  way,  you  can  add  devices  using  higher  numbers,  and  not  affect  the 
;  devices  already  defined.  Put  your  favorite  output  formats  at  the  top. 

device_number  =  [0,  9,  2,  X,  3,  4,  10,  5,  6,  7,  8] 


temp  = 

t  =  0 

print, 

print, 

print, 

print, 

print, 

print, 

print, 

print, 


1. )  VT340  Screen' 

2. )  PostScript  Printer' 

3 .  )  X  Window' 

4 .  )  QMS  Printer' 

5. )  Tek4014  Screen' 

6. )  VT340  File' 

7. )  PostScript  File' 
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print,  '  8.1  QMS  File1 

print,  '  9.)  Tek4 014  File' 

print,  '  10.1  HPGL  File' 

if  keyword_set  ( polar __opt ion)  then  begin 

print,  '  11.1  Plot  (pre-printed  paper  format  on  TALAR I S 1 ) ' 

endif 

print,  '  ' 

read,  'Choose  format  for  graphics  output:  temp 

if  temp  ne  ''  then  t  =  fix (temp) 
if  t  ge  1  and  t  le  9  then  begin 
device  =  device_number ( t  -  1) 
device_name  =  device_names (device) 
endif 
return 
END 


pro  3etup_device,  limited 


;  This  routine  should  be  customized  so  that  files  are  sent  to  the  proper 
;  printers.  Feel  free  to  add  new  device  codes.  You  should  also  take  a 
;  look  at  reset_device .pro  if  you  are  going  to  make  modifications. 

common  devcom,  device,  device_name,  vx,  wy,  orient 

CASE  device  OF 

;  VT340  screen 
0:  BEGIN 

set_plOt,  'REGIS' 

END 

QMS  printer 
1:  BEGIN 

print,  'Creating  PLOT. QMS  and  sending  it  to  TA1ARIS1 . . . ' 
set_plot,  'QMS' 
if  orient  eq  l  then  begin 

device,  f ilenamer 'plot .qms ' ,  /landscape,  xoffset=.25,  $ 
yof f set= . 25 , xsize=12 . 0 , ysize=12 .0 , /inches 
endif  else  begin 

device,  filenames 'plot .qms ' ,  /portrait,  xoffset=.25,  $ 
yof f sets ,25 ,xsize=12 .0,ysize=12 .0, /inches 

endelse 

if  not  keyvord_set ( limited)  then  begin 
ip. thick  =  1.1 
<x. thick  =  1.1 
iy. thick  s  1.1 
ty. margin  =  (4,31 

tx  =  (wx  *  2.54  *  !d.x_px_cm}  /  *d.x_vsize 
ty  =  (wy  *  2.54  *  -d.y_px_cm)  /  !d.y_vsize 
if  orient  eq  1  then  tx  =  tx  *  1.012  else  ty  =  ty  *  1.012 
set_vievport,  (l-txl/2,  l-(l-tx>/2.  U-ty>/2,  l-U-tyl/2 
endif 

END 

;  X  Window 
2:  BEGIN 

8et_plot,  'X' 

END 

;  Tek4014  screen 
3:  BEGIN 

set^plot,  'TEK' 

END 

;  VT340  graphics  file 
4:  BEGIN 
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print,  'Creating  PLOT. REG...' 
set_plot,  'REGIS' 
get_lun,  lu 
openv,  lu,  'PLOT. REG' 
device,  plot_to=lu 

END 

QMS  graphics  file 
S:  BEGIN 

print,  'Creating  PLOT. QMS...' 

set_plot,  'QMS' 

if  orient  eq  1  then  begin 

device,  t ilename= 'plot .qms ' ,  /landscape,  xoffset=.25,  $ 
yof f set- .25, xsize=l2 .0,y sizes 12 .0, /inches 
endif  else  begin 

device,  t ilename^' plot .qms' ,  /portrait,  xoffset=.25,  $ 
yof f set= .2S,xsi2e=l2 . 0,ysize=12 .0, /inches 

endelse 

if  not  keyvord_set ( limited)  then  begin 
ip. thick  =  1.1 
!x. thick  =  1.1 
)y. thick  =  l.i 
:y. margin  =  (4,31 

tx  =  (wx  *  2.54  *  id.x_px_cm)  /  !d.x_vsize 
ty  =  I wy  *  2.54  *  !d.y_px_cm)  /  id.y^vsize 
if  orient  eq  1  then  tx  =  tx  *  1.012  else  ty  =  ty  *  1.012 
set_viewport, (l-tx>/2,  l-(l-tx)/2,  (l-ty}/2,  l-(l-ty)/2 
endif 

END 

Tek4014  graphics  file 
6;  BEGIN 

print,  'Creating  PLOT-TEK...' 
set_plot,  'TEK' 
get_lun,  lu 
openw,  lu,  'PLOT. TEK' 
device,  plot_to=lu 

END 

HPGL  graphics  file 
7:  BEGIN 

print,  'Creating  PLOT.hpg...' 
set_plot,  'HP' 

device,  filenames 'plot .hpg' ,  /landscape,  /inches 

END 

device  8  is  prints  a  special  polar  plot  file  to  the  QMS  printer.  This 
plot  can  be  over lay ed  on  top  of  special  polar  graph  paper  that  we  use 
for  measurements  so  we  can  compare  modeled  data  with  actual  measurements. 

PostScript  printer 
9:  BEGIN 

print,  'Creating  PLOT. PS  and  sending  it  to  POST$M0DEL_LNO5 . . . ' 

set_plot,  'PS' 

if  orient  eq  1  then  begin 

device,  filenames 'plot .ps' ,  /landscape,  xof£set=0.0,  $ 
yof fset=ll .0 ,  xsize=ll ,0,ysize=8.5, /inches 
endif  else  begin 

device,  filenames 'plot .ps' ,  /portrait,  xoffsetsO.O,  $ 
yof fset=0 .0,xsize=8 .5,ysize=ll .0, /inches 

endelse 

if  not  keyvord_set ( limited)  then  begin 
;p. thick  =  1.1 
!x. thick  =  1.1 
•y. thick  =  1.1 
!y. margin  =  [4,3] 

tx  =  (wx  *  2.54  *  !d.x_jpx_cm)  /  id.x^vsize 
ty  =  (vy  *  2.54  *  !  d  .y _px_cm)  /  !d.y_vsize 
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set .viewport , ( 1-cx) /2 ,  l-[i-tx)/2.  (i-ty)/2.  l-(l-Cyi/2 

print,  (1-cx)  /2.  l-U-tx)/2.  (l-cy)/2.  l-(l-ty)/2 
hak 
endif 

END 

;  PostScript  graphics  file 
10:  BEGIN 

print,  'Creating  PLOT. PS...' 

set_plot,  'PS' 

if  orient  eq  1  then  begin 

device,  f ilename= 'plot .ps ' ,  /landscape,  xoffset=0.0,  $ 
yof fset=0 . 0,xsize=ll .0,ysize=8 .5, /inches 
endif  else  begin 

device,  filenames ‘plot .ps ' ,  /portrait,  xoffset=0.0,  $ 
yof f set= 0 . 0,xsize=8 . 5,ysize=11.0, /inches 

endelse 

if  not  keyword_set ( limited)  then  begin 
ip.  thick.  =  1.1 
tx. thick  =  1.1 
iy. thick  =  1.1 
iy. margin  =  (4,31 

tx  =  <wx  *  2.54  *  !d.x_px_cm)  /  id.x_vsize 
ty  =  (wy  *  2.54  *  id.y_px_cm)  /  !d.y_vsize 
set_viewport , ( 1-tx) /2 ,  l-(l-tx)/2,  (l-ty)/2,  l-(l-ty)/2 
endif 

END 

ELSE:  ;  do  nothing 

ENDCASE 
return 
end 


pro  reset_device 


;  This  routine  resets  variables  and/or  closes  files,  setting  the  plot 
;  device  to  the  default  device  (REGIS  1 .  Feel  tree  to  modify  this  routine 
,*  in  order  to  use  your  particular  graphics  devices  properly.  Be  sure  to 
;  look  at  setup_device.pro  also. 

common  devcom,  device,  device_name,  wx,  wy,  orient 

CASE  device  OF 

;  VT340  screen 
0: BEGIN 

set_cursor,  23,  20 
hak,  /mesg 
9et_pIot,  'REGIS' 

END 

;  QMS  printer 
1 : BEGIN 

ip. thick  =  0 
*x. thick  =  0 
iy. thick  =  0 
iy. margin  *  (4, 2] 
device,  /close 
setjplot,  'REGIS* 
spawn,  'cprint  plot.qms' 
print,  'Done.' 

END 

;  X  window 
2: BEG IN 

set^plot,  'REGIS' 
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END 

;  Tek4014  screen 
3 : BEGIN 

set_cursor,  23,  20 
hak,  /mesg 
set_plot,  'REGIS' 

END 

;  VT340  graphics  file 
4 ; BEGIN 

print,  'Done.' 
device,  /close,  plot_to=0 
close,  lu 
set_plot,  'REGIS' 

END 

;  QMS  graphics  file 
5: BEGIN 

i p . thick  =  0 
!x. thick  =  o 
!y. thick  =  o 
!y. margin  =  (4,2] 
print,  'Done.' 
device,  /close 
set_plot,  'REGIS' 

END 

;  Tek4014  graphics  file 

6: BEGIN 

print,  'Done.' 
device,  /close,  plot_to=0 
close,  lu 
set_plot ,  'REGIS' 

END 

;  HPGL  graphics  file 

7: BEGIN 

print,  'Done.' 
device,  /close 
set_plot,  'REGIS' 

END 

9: BEGIN 

!p. thick  =  o 
:x. thick  s  0 
!y. thick  =  o 
!y. margin  =  (4,2] 
device,  /close 
set_plot,  'REGIS' 
spawn,  'post  plot .ps ' 
print,  'Done.' 

END 

10 : BEGIN 

!p. thick  =  0 
!x. thick  =  0 
iy. thick  =  0 
sy. margin  =  (4,2] 
device,  /close 
set_plot,  'REGIS' 
print,  'Done.' 

END 

ELSE : 

ENDCASE 

return 

end 
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Appendix  B 

SUBROUTINE  DESCRIPTIONS 


This  appendix  contains  a  brief  description  of  each  subroutine  and  a  table  (Table  B 1 )  showing  the 
interdependence  of  the  various  subroutines.  The  hierarchy  of  subroutines  is  from  left  to  right  in  the 
table. 


append_file:  Appends  the  file  and  data  names  of  a  new  data  object  to  the  structures 

maintaining  this  information. 


append_p_info: 

cfit_interface: 

color_select: 

gauss2: 

genload: 


graphics; 

histgram2: 

histo _plot: 
laminate: 


legend: 


Appends  the  default  plotting  parameters  and  data  pairs  to  the  structures 
maintaining  this  information  for  the  standard  and  polar  single  plot  routines. 

Allows  the  user  to  select  one  of  several  different  curve  fitting  options  and 
then  performs  the  fit  of  the  data  to  the  specified  model. 

Allows  the  user  to  select  a  color  for  a  plot  line. 

Performs  gaussian  elimination  with  maximal  pivoting  the  rational 
curvefitting  routine. 

A  general  purpose  function  to  load  data  files  into  a  variable.  For 
unformatted  files,  it  asks  the  user  about  data  types  and  number  of  columns. 
For  formatted  files,  the  number  of  columns  and  lines  of  data  are  determined 
automatically. 

The  main  PV-WAVE  program.  Most  of  the  higher  level  menus  arc  here. 

Sets  up  and  plots  statistical  information  about  a  data  set  (histogram,  mean, 
median,  mode,  standard  deviation,  range,  min,  and  max). 

Main  menu  controlling  the  statistical  plots. 

Appends  a  vector  to  an  array.  If  the  vector  is  smaller  than  the  array,  then  it 
is  padded  with  zero's  to  match  the  length  of  the  array.  If  the  array  is 
smaller,  then  the  array  is  extended  in  size  by  padding  with  zeros  to  match 
the  length  of  the  vector. 

The  legend  menu  is  here.  It  allows  the  user  to  define  a  legend  for  standard 
and  polar  overplots.  The  user  has  control  over  the  position  of  the  legend  on 
the  graph,  text  for  each  plot  entry,  line  type,  and  color. 
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line_descnpt: 

line_menu. 

line_select: 
loaded_.fi  les: 

loader 


Lists  the  line  styles,  colors,  and  legends  defined  for  all  marked  data  pairs  to 
be  used  in  an  overplot  operation. 

Allows  the  user  to  select  a  line  r  yle,  color,  oi  legend  for  a  marked  data  pair 
in  an  overplot.  This  routine  calls  line_descript  to  list  the  default  values. 

Allo  ws  the  user  to  select  a  line  style  for  a  plot  line  in  the  histogram  section. 

Prints  a  list  of  the  files  loaded  into  the  program.  Each  entry  is  separated  by 
spaces 

Asks  the  user  for  a  file  name  to  load.  It  calls  genload  and  can  print  out  a  list 
of  data  files  in  the  current  directory. 
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manip_help. 

manipulate: 

mode: 

mouse_utilities: 

name_check: 

output_select: 

polar_oplot: 

polar_plot: 

polar_pnnt: 

printer_parms: 

ratfit: 

ratsetup: 

save_data: 

selea_hst: 

select_manip: 

select_multi: 

select_plot: 

show_manip_info: 

show_selected: 


Lists  the  operations  and  functions  supported  by  the  manipulate  operation. 

Allows  data  to  be  created  by  applying  mathematical  operations  to  existing 
data. 

Computes  the  statistical  mode  of  a  set  of  data 

Provides  a  series  of  options  within  the  standard_plot  menu  that  require  a 
mouse.  These  options  include  finding  the  min/max  of  a  data  range,  finding 
the  closest  point,  and  zoom. 

Prompts  for  a  data  name  for  newly  created  data  and  ensures  that  the  name 
provided  is  unique. 

Allows  the  user  to  select  a  destination  device  for  a  plot. 

Polar  overplot  menu  and  routines 

Plots  one  dataset  assuming  that  the  data  is  in  polar  coordinates. 

Specialized  polar  print  routine  that  scales  the  plot  so  that  it  can  be  overlaid 
on  standard  NRL  polar  graph  paper. 

Allows  the  user  to  define  the  width  and  height  in  inches  and  whether  the 
printout  will  be  in  landscape  or  portrait  mode.  Currently  works  only  with 
QMS  printers  but  should  work  with  postscript  printers  also. 

Performs  rational  fraction  curvefitting. 

Allows  the  user  to  select  the  coefficients  for  fitting  in  the  rational  fraction 
curvefitting. 

Writes  a  newly  created  data  object  to  a  file. 

Allows  the  user  to  mark  any  number  of  single  data  sets  and  to  select  one  of 
these  data  sets  for  display  in  a  histogram  (statistical  plot). 

Allows  the  user  to  select  one  or  two  data  columns  to  use  for  data 
manipulation.  The  user  may  also  select  a  previous  selection  for  re-use. 

Similar  to  select_plot  It  allows  the  user  to  select  multiple  sets  of  data  that 
will  be  plotted  together  by  either  the  overplot  or  multiplot  functions. 

Allows  the  user  to  create  any  number  of  paired  data  columns  and  then  select 
one  of  these  data  sets  for  an  upcoming  operation. 

Lists  the  data  that  has  been  previously  selected  for  data  manipulation. 

Lists  the  data  column  pairs  that  have  previously  been  marked  as 
ordinate/abscissa  or  angle/range  pairs. 
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standard_mplot: 

standard_oplot: 

standard_plot: 

xypick: 


Standard  multiplot  menu  allows  the  user  to  display  more  than  one 
independent  plot  on  one  or  more  separate  pages/screens. 

Standard  overplot  menu  plots  multiple  data  sets  on  one  graph. 

Standard  plot  menu  only  plots  one  data  set. 

Allows  the  user  to  select  a  column  of  data  from  one  of  the  data  sets 
currently  loaded. 
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Table  B1  -  Subroutine  Hierarchy 


graphics 


loaded_filea 
loader  - 


-4  genload 


append_file 

c£it_interface 


►  sare_data 

•  name_cheek 

►  selectplot 

manipulate  - 
xrpick  - 


r 


ratf it 


name  cheek 


t: 


4  gauss 2 


ratsetup 


-4  show_selected 
■4  *anip_help 
■4  loaded  files 


x: 


polarplot - r~4  polarprint 

— 4  printerparas 
— ►  printerpetup 
— 4  printerpeset 

►  standardplot^J^  »  «ouse_utilities 

►  selectpulti  4  shawpelected 

’  pol^r^oplot  -  |  »  polarprint 


z 


lineaenu 


U 


4  line_descript 
showpe lected 


standardpplot  _n 

selectjist  - 

his  top  lot 


printerpj 
setup_derice 
u  reset_device 
4  show  hstselect 


4  histgraa2  - 

4  colorpelect 
4  linepelect 
4  setup_deTi.ee 
►  reset_derice 


-4  mode 


— *  laminate 

1—4  standardpplot 
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